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Abstract 



We propose an iterated Azema-Yor type embedding in the spirit of 



Azema and Yor, 19791 and Brown et al., 2001a for any given finite num- 



ber n G N of probability measures //!,...,//„ which are in convex order and 
satisfy an additional technical assumption. In particular, our construction 



reproduces the stopping boundaries obtained in Madan and Yor, 2002 



and Brown et al., 2001a 



We demonstrate with a counterexample that 
our technical assumption is necessary and propose an extended embed- 
ding for n = 3. As a by-product of our analysis we compute the law 
of the maximum at all the stopping times. This is used in our parallel 



work iHenry-Labordere et al., 2013 to establish extremal properties of our 



embedding and develop robust pricing and hedging of Lookback options 
given prices of call options with multiple intermediate maturities. 

1 Introduction 

The Skorokhod embedding problem (SEP) consists in finding a stopping time r 
such that a given stochastic process (Xj) stopped at r has a given distribution 
fi. For the solution to be useful (and non-trivial) one further requires r to be 
minimal (cf. [Obloj, 2004 Sec. 8]) which, when X is a continuous local martingale, 
is equivalent to {Xt^r '■ t ^ 0) being a uniformly integrable martingale. The 
problem, dating back to the original work in Skorokhod, 1965 , has been an 
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active field of researcli for nearly 50 years. New solutions often either considered 
new classes of processes X or focused on finding stopping times r with additional 
optimal properties. In particular, the latter is true for works on SEP motivated 
by mathematical finance. 

The classical approach in quantitative finance postulates a model - a stochas- 
tic process on a filtered probability space - and uses the model together with 
no-arbitrage assumption to deduce the fair price today of an option with pay- 
off 0{X)t, where 0{X)t is simply an F^-measurable random variable. The 
no-arbitrage assumption implies, up to an equivalent change of measure and dis- 
counting, that our stochastic process is a martingale and then the price is given as 
E [0(X)t]. However in practice many of such options have known prices quoted 
in the market. It is therefore natural to pose an inverse problem: given a set of 
options with market prices what is the range of possible prices for 0{X)j' among 
all consistent models? This logic is often referred to as robust, or model-free, ap- 
proach. If we assume that market quoted prices for options with payoff (Xt — K)"^ 
(call options) are given for all K then this is equivalent to fixing the marginal 
law Xt ~ fi- If we further assume that X has continuous paths then consistent 
models correspond to all continuous martingales with a given marginal at time 
T. These can naturally be represented as time-changes of Brownian motion B 
and, as long as 0{X)t is invariant under time change, the robust pricing question 
above boils down to finding the range of E[0(i?)T-] among all solutions to the 
Skorokhod embedding problem. 

The above link between SEP and robust pricing and hedging was pioneered 



by Hobson, 1998 who considered Lookback options. Barrier options were sub- 
sequently dealt with by Brown et al., 2001b . More recently, |Cox and Qbloj, 



rier options. 



2011b Cox and Obloj, 2011a considered the case of double touch/no-touch bar- 



Hobson and Neuberger, 2012 looked at forward starting straddles 

We 



and analysis for variance options was undertaken in Cox and Wang, 2013 



refer to [Hobson, 2010 and Obloj, 2010 for an exposition of the main ideas and 
more references. 

All the above works assumed that prices of call options with only one maturity 
are given, or equivalently that one marginal at the final time T is fixecjj Mean- 
while, in practice, often prices of call options with several intermediate maturities 
are also available. Mathematically, this translates into optimising E [0(-B)t-„] over 
all solutions ri < • • • < r„ to a multi-marginal Skorokhod embedding problem: 



/ij, 1 < i < ra, and {B, 



t^TnJt^Q 



is a uniformly integrable martingale, where 



Br 

the convex order. As shown by Strassen, 1965 , the convex ordering of measures 



, /i„ is a given sequence of centred probability measures non-decreasing in 



is both necessary and sufficient for a solution to exist. 



^With the exception of [Hobson and Neuberger, 2012 where a second marginal corresponding 
to the starting date of the option is also fixed leading to a one-marginal SEP but with a non- 
trivial starting law. 



This problem turns out to be significantly more complex that the special case 
n = 1. In fact, we are not aware of any previous works which would solve it 



for non-trivial examples of O. The only exception is Brown et al., 2001a where 
the authors constructed a generalisation of the Azema-Yor embedding to two 
marginals and proved its optimal properties thus solving the problem for n = 2 
and 0{B)r = 0(supj^^i?j) for an increasing function (p. In this paper we extend 
this further and treat the case of arbitrary n e N. We prove the embedding 
property for a general sequence of measures which satisfy an additional technical 
assumption. 

We started our quest for a general n-marginal optimal embedding by comput- 
ing the value function sup E [^(supj^^^ 5^)] among all solutions to the n-marginal 
SEP. This was achieved using stochastic control methods, developed first forn = 1 
in Galichon et al., 2013 , and is reported in a companion paper Henry-Labordere 
et al., 2013 . Knowing the value function we could start guessing the form of the 
optimiser and this led to the present paper. Consequently the optimal properties 
of our embedding, namely that it indeed achieves the value function in ques- 



tion, are shown in Henry-Labordere et al., 2013 . In fact we give two proofs in 



that paper, one via stochastic control methods and another one by constructing 
appropriate pathwise inequalities. 

The reason why the n-marginal problem is much more involved than the one- 
marginal one is that while it might be easy to iterate a one-marginal solution, 
there is no guarantee that such a simple iteration of optimal embeddings would 
be globally optimal. Indeed, this is usually not the case. Consider the Azema- 



Yor embedding |Azema and Yor, 1979 which consists of a first exit time for the 
joint process {Bt,Bt)t^Q, where Bf = sup^^^-Bs- More precisely, their solution 
T^^ = inf |i ^ : -Bt < ^^{Bt)] leads to a functional relation B^-ay = ^^{B^ay). 
This then translates into the extremal property that the distribution of B^ay is 
maximized in stochastic order amongst all solutions to SEP for fi, i.e. for all y, 

P [B^AY ^y]= sup {P [Bp ^y]:p s.t. B^ ^ /i, (5^,,) is UI } . 

It is not hard to generalise the Azema-Yor embedding to a non-tirivial start- 
ing law, see Obloj, 2004[ Sec. 5]. Consequently we can find T]i such that Tj = 



inf |t ^ rj_i : Bt ^ rji{sup^.__^^g^^ Bg)^ solve the n-marginal SEP. However this 
construction will maximise stochastically the distributions of sup^^^^^^^^ Bt but 
not of the global maximum B-r^ . 

Instead, we propose an "iteration" of the Azema-Yor embedding in the follow- 
ing sense. Each Tj+i is still a first exit for {Bt,Bt)t^T-. which is designed in such 
a way as to obtain a "strong relation" between -Bri+i and -Bri+u ideally a func- 
tional relation. Under our technical assumption about the measures /ii, . . . ,/i„. 



Assumption |A| we describe this relation in detail in Lemma 3A which is key 
for the pathwise proof of the optimal properties of our embedding in Section 4 
in [Henry-Labordere et al., 2013 . 



For n = 2 our work recovers the results in Brown et al., 2001b . We also 



recover the trivial case Tj = r^^ which happens when ^^. < ^/^,+i, we refer to 



Madan and Yor, 2002 who in particular then investigate properties of the arising 
time-changed process. However, as a counterexample shows, our construction 
does not work for any laws /xi, . . . , /i„ which are in convex order. A more involved 
analysis is required to obtain an explicit embedding when our Assumption [Ajfails. 
We only sketch the appropriate arguments for the case n = 3 as this case does 
not seem to be of any practical relevance. 

The remainder of the article is organized as follows. In Section |2] we explain 
the main quantities for the embedding and state the main result. We also present 
the restriction on the measures /xi, . . . , /i„ which we require for our construction 
to work (Assumption IaI) . In Section [3] we prove the main result and Section |4 
provides a discussion of extensions together with comments on Assumption |Aj 
The proof of an important but technical lemma is relegated to the Appendix. 



2 Main Result 

Let {Q, J-", F, P) be a filtered probability space satisfying the usual hypothesis 
and B an F-continuous local martingale, Bq = 0, (-B)oo = oo a.s. and B has 
no intervals of constancy a.s. We denote Bt := sup^^^-Bt. We are primarily 
interested in the case when 5 is a standard Brownian motion and it is convenient 
to keep this example in mind, hence the notation. We allow for more generality 
as this introduces no changes to the proof of statements. 

2.1 Definitions 

The following definition will be crucial in the remainder of the article. We define 
the stopping boundaries ^i, . . . ,^„ for our iterated Azema-Yor type embedding 
together with quantities Ki, . . . , Kn which will be later linked to the law of the 
maximum at subsequent stopping times. 

Definition 2.1. Fix n g N. For convenience we set 

co^O, i^o^O, i^{y) = -^. (2.1) 

For C ^ H^ o-nd i = 1, . . . ,n we write 

q(C):= f ix-Cffii{dx). (2.2) 



Let y ^ and assume that for i = 1, . . . ,n — 1 the quantities ^i,Ki, li and ji 
are already defined. Then we define 

in{-;y) : {-co,y] ^ {0,l,...,n-l}, 

«n(C; y) ■= max {A; e {0, 1, . . . , n - 1} : ^k{y) < C} , 



and 



iniy) := sup ^ arg inf \ -^^— 



g»n(C;j/)(C) 



K. 



*n(C;3/) 



(y) 



With 

we set 
Kn{y) : 



]n{y) ■=in{^niy);y) 



(2.4) 



(2.5) 



^^niUy)) - [cMy)iUy)) -{y- Uy))KMy)iy)] ■ (2.6) 



y - in{y) 

Definition 2.2 (Embedding). Set r = and /or z = 1, . . . ,-«. define 
inf {t ;^ r,_i : 5, ^ ^^(5*)} if B,^_^ > UK-,)^ 



Ti 



Ti-l 



else. 



(2.7a) 
(2.7b) 



Remark 2.3 (Properties of in) ■ Recalling the definition of in, cf. (2.3 ), we observe 
for later use that for y ^ 



tn{-',y) is left-continuous and has at most n ~ 1 jumps 



and for x e 



in{x] ■) is right-continuous and has at most n — 1 jumps. 



(2i 



(2.9) 



Figure [23] illustrates a set of possible stopping boundaries in the case of n = 3. 
If Assumption |X] is in place, see Section 2^, we will show that the stopping 
boundaries are continuous (except possibly for i = 1) and non-decreasing, cf. 
Section 12. 5[ 

The n**^ stopping boundary ^„ is obtained from an optimization problem which 
features ^i, . . . , ^n-i and Ki, . . . , Kn-i- Kn{y) is the value of the objective func- 
tion at the optimal value C,n{y)- Note that all previously defined stopping bound- 
aries ^1, . . . , Cn-i and the quantities Ki, . . . , Kn-i remain unchanged. 

Denote the right endpoint of the support of the measure /ij by 



r^. := inf {x : /i^ {{x, oo)) = 0} 
and the barycentre function of /ij by 



[U 



'^ ' -i([x,oo)) ^<^>J I^^^mJ' 



(2.10) 



(2.11) 




Figure 2.1: We illustrate possible stopping boundaries ^1,^2,^3- The horizontal 
lines represent a sample path of the process [Bt, Bt) where the x-axis is the value 
of B and the y-axis the value of B. Each horizontal segment is an excursion of B 
away from its maximum B. The first stopping time Ti is found when the process 
first hits ^1. Since ^i{BrJ > ^2(-Bti) the process continues and targets ^2- The 
stopping time T2 is found when the process first hits ^2- Since ^2(-Br2) ^ ^si^r^) 
we get Ts = T2. For the y we fixed we have i3{xi; y) = 0, Z3(x2, y) = 1, i-iix^; y) = 2. 



As shown by Brown et al., 2001b , the right-continuous inverse of 6j, denoted 
by hj^, can be represented as 



h-\y) 



sup < arg mt 

<:<y y-C 



(2.12) 



It is clear and has been studied in more detail by Madan and Yor, 2002 that 



if the sequence of barycentre functions is increasing in i, then the intermediate 
law constraints do not have an impact on the corresponding iterated Azema-Yor 
embedding. However, in general the barycentre functions will not be increasing 
in i, cf. [Brown et al., 2001a |, and hence will affect the embedding. We think 
of jn{y) as the index of the last law /ij,i < n, which represents, locally at level 
of maximum y, a binding constraint for the embedding. As compared to the 



optimization from which 6„ is obtained, cf. (2.12), the optimization from which 



^n is obtained, cf. (2.4), has a penalty term. 



2.2 Restrictions on Measures 



Throughout the article we will denote the left- and right-limit of a function / at 
X (if it exists) by /(x— ) and /(x+), respectively. 

For ra e N and y ^ define inductively the mappings 



c"(-,i/) :(-(», 1/]^M 



X ^-^ c 



\x,y) := Cnix) - [cr,^{x-y)ix) - (y - x)K^^(x.y){y)] 



(2.13) 



It follows that the minimization problem in (2.4) is equivalent to the following 
minimization problem, 



U{y) e argmm — , 



(2.14) 



where we observe that 



C'iCy) 



y-C 



<:=y 



lim^^ (2.15) 

av y-C 

-C'niy-) + C[^(y;y){y-) + K^„{y;y){y) if Cn{y) = C^„{y,y){y), 

+00 else, 



and note that existence in (2.14) can be deduced from the - a priori - piecewise 
continuity of c"'(-, y), the fact that c" ^ and the property that C ^^ ^'^' — ^^^ 
for ( sufficiently small, which is a non-increasing function. 
For later use observe 



y-( 



y-c. 



minfc^ {y) ^ ^„(y) ^ y 



(2.16) 



which follows from the definition of ^„, cf. (2.4), and where b^ denotes the 



right-continuous inverse of the barycentre function 6j, cf. (2.12). 



Assumption A (Restriction on Measures). Recall definitions in (2.1)-(2.2), 



(2.13) and (2.10). We impose the following restrictions on /ii, . . . ,/i„.- 



(i) 5 |a;|/ij(dx) < oo with '\xfii{dx) = and Cj_i < q for all 1 ^ i ^ n, 
(ii) for all 2 ^ i ^ n and all ^ y ^ r^- the mapping 

cUC,y) 
(—00,2/] ^- K u 1+°'^}) C *~^ ^ has a unique minimizer (* (2.17) 



y-C 



and 



Ci{y) > c^^^y.y){y) whenever C < y- 



(2.18) 



Remark 2.4 (Assumption^ . The condition that the call prices are non-decreasing 
in maturity 



Ci < c 



J + l5 



1 n 



(2.19) 



can be rephrased by saying that /ii , . . . , /i^ are non-decreasing in the convex 
order. Condition (i) in the Assumption is the necessary and sufficient condition 
for a uniformly integrable martingale with these marginals to exist, as shown by 
e.g. 



Strassen, 1965 



Note that if (2.19) holds with strict inequality then (2.18) is automatically 
satisfied. 



Remark 2.5 (Discontinuity of ^i). Note that Assumption jAKii) does not require 
that the mapping 



c 



y-C y-C 



(2.20) 



has a unique minimizer. It may happen that there is an interval of minimizers 
and then ^i is discontinuous at such y. 



2.3 Main Result 

Our main result shows how to iteratively define an embedding of (/xi, . . . , /i„) in 
the spirit of Azema and Yor, 1979 and Brown et al., 2001a if Assumption [A| is 
in place. 

Theorem 2.6 (Main Result). Let n e N and fii, . . . , fin be probability measures 
on M. which satisfy Assumption^ in Section \27^ Recall Definitions 2A_ and \2.^ 
Then Ti < oo, Br^ ~ yUj for alii = 1, . . . ,n and {Br„^t)t>o '^^ ^ uniformly integrable 
martingale. 

In addition, we have for y ^ and i = 1, . . . ,n, 



F[B^^^y]=K,{y) 



(2.21) 



where Ki is defined in (2.6) 



Remark 2.7 (Inductive Nature). It is important to observe that ^j and therefore 
also Tj, only depend on /xi, . . . ,/ij. This gives an iterative structure allowing to 
"add one marginal at a time" and enables us to naturally prove the Theorem by 
induction on n. 

Remark 2.8 (Minimality). Since all Tj are such that {Bt^T--)^^^ is a uniformly 
integrable martingale it follows from Monroe, 1972 that all Tj are minimal. 



2.4 Examples 



Examples |2.9|and|2.10 , respectively, show that we recover the stopping boundaries 
obtained in Madan and Yor, 2002 and Brown et ah, 2001a , respectively. In 



particular the case n = 1 corresponds to the Azema-Yor solution Azema and 



Yor, 1979 



Example 2.9 ( Madan and Yor, 2002|). Recall the definition of the bary centre 
function hi from (2.11). 'Madan and Yor, 20021 consider the "increasing mean 
residial value" case, i.e. 



6i < 62 < • • • ^ &n. 



(2.22) 



We will now show that our main result reproduces their result if Assumption 
pjzfi in place. In fact, as can he seen below, our definitions of ^i and Ki, cf. (2.4) 
and (2.6), respectively, reproduce the correct stopping boundaries in the general 
case, showing that Assumption [^ is not necessary, cf. also Section |^ More 
precisely, we have 



ii 



(,-1 



K.{y) - ^^^KM\ -■ ^^fWvMl 



(2.23) 



y~K\y) 
where b^^ denotes the right- continuous inverse of hi and nf^ is the Hardy-Littlewood 



transform of Hi, cf. ^Carraro et al., 20121 



Clearly, the claim is true for i = 1. Let us assume that the claim holds for all 



2 < n — 1. Now, the optimization problem for ^„ in (2.4) becomes 



e„(y)6argmin|^-l^^^,-.^(^)| 



Cn-l(C) C„_i(&„_i(y)) 

y^C y-b-\{y) 



arg mm < mm 



Cn(C) 



mm 



c^y c^b-\{y) y-C c^Cife) \y-C 



Cn(C) 



Cn-l(C) _ Cn-l{h^\{y)) 

y-C y~b-\{y) 



It is clear that the first minimum is Ai = .-i/^! since b^^{y) < 6„\(y). 
As for the second minimum, we set 

F(C) = MO _ (c_n-i{0 ^ c4b-\{y)) 

^'' y-C [y-C y-b;Uy) 

and we see by direct calculation that 

{y - C?F'{C) = (fen(C) - 2/)^n ([C, »)) - (&n-l(C) y) ^^n-l ([C, »)) 

fy-, bn{0-y bn-i{C)-y 

= Cn(C)rT7\ 7 - C„_i(C): 



bniO C 



bn^liO-C 



where derivatives exist for almost all ( ^ M.. By (2.22), we conclude therefore 

y 



iy-CYF'iO^icniO-Cn-iiO) 



M-l 



(C) 



&n-l(C)-C 



^0, 



where the last inequality follows from the non-decrease of the Hi 's in the convex 
order. Hence F is non- decreasing, and it follows that it attains its minimum 



at the left boundary, i.e. A2 

Consequently, niin{y4i, ^42} = Ai and (2.23) follows 



c,+i{b-'{y)) 






c,+i(6-i)(y) 

y-b-\y) 



Example 2.10 ( [Brown et al., 2001a|). In the case of n = 2 our definition of ^ 



and ^2 clearly recovers the stopping boundaries in the main result of ^ Brown et al. 



2001 ah However, our embedding is not as general as their embedding because we 



enforce Assumption\^ see also the discussion in Section^ 
Example 2.11 (Locally no Constraints). In general we have 

However, if 

Uy) = K\y) 



for some y ^ then it follows from Theorem \2.6\ that 

Cn{b-\y)) 



Kn{y) 



y-K^{y) 



^^T{[yMl 



(2.24) 



(2.25) 



(2.26) 



i.e. locally at level of maximum y the intermediate law constraints have no impact 
on the distribution of the terminal maximum as compared with the Azema- Yor 
embedding. 

2.5 Properties of ^„ and Kn 

Under Assumption |X] we establish the continuity of C,n for n ^ 2, cf. Lemma 
2.12 and prove monotonicity of ^„ for n ^ 1, cf. Lemma 2.13 In Lemma 2.14 
we derive an ODE for Kn which will be later used to identify the distribution of 



the maximum of the embedding from Definition 2.2 



Let Til < 712- Recalling Remark 2.7 it follows that the embedding of the 
first Hi marginals in the n2-marginals embedding problem coincides with the 
rii-marginals embedding problem. Hence it is natural to prove the Lemma by 
induction over the number of marginals n. 

Lemma 2.12 (Continuity of ^„). Let n^ 2 and let Assumptions^ hold. Set 

A := {(x,i/) eM X M+ : x<i/}. (2.27) 



10 



Then the mappings 



c": A 



y ^ ^n{y) 



(2.28) 
(2.29) 



are continuous. 



Proof. We prove the claim by induction over n. Let us start with the induction 
basis n = 1,2. Continuity of c^ is the same as continuity of Ci and continuity of 



c^ is proven in 



Brown et aL, 2001a , cf. Lemma 3.5 therein. As for continuity of 



^2 we note that our Assumption Ajii) precisely rules out discontinuities of ^2 as 



seen in Section 3.5 of [Brown et al., 2001a . By induction hypothesis we assume 
continuity of c^, . . . , c""^ and ^2, • • • , Cn-i- 

The only possibility that a discontinuity of c" can occur is when the index tn 
changes. This only happens at {x,y) = {C,k{y),y) for some k < n, or, in the case 
that y is a. discontinuity of ^1, at {x,y) where x e [^i{y^),C,i{y+)]. We prove 
continuity at {x,y). 

Consider first the following cases: 



or, 



if X = ^kiy) 

ifxe[ei(y-),6(y+)] 



then X ^ Cj{y) for all j y^ k, j < n, (2.30) 
then X ¥= C,j{y) for all j ¥= I, j < n. (2.31) 



Note that in case (|2.31[) we have from Remark |2. 5 



KM 



CiiX) 



X 



for all XE [ii{y-),ii{y+)]. 



(2.32) 



We will call a point {x,y) to be "to the right of ^a:" if ikiy) < x and "to the 



left of ^fc" ii^k{y) ^ X. From (2.30) and (2.31) it follows that there exists an e > 
such that each point (x, y) in the e-neighbourhood of (x, y) is either to the left 
or to the right of ^k and there are no other boundaries in this e-neighbourhood, 
in particular 



k 



'"nyXTi Vr ) y 



J '^n[Xl,yi) \n{xr;yr)\Xr)yr)i 



(2.33) 



where (x^, yr) is in the e-neighbourhood of (x, y) and to the right of $,k and (x^, yi) 
is in the e-neighbourhood of (x, y) and to the left of ^k- 



11 



li X < y, we have by induction hypothesis 

C^{Xr,yr) = Cn{Xr) - {Cfc(Xr) " {Vr - Xr)Kk{yr)} 

Cn{x) - {cfc(x) ~{y- x)Kk{y)) 
y-x 



(2.34) 



{xr,yr)-^{x,y) 
> 

from the right 



c„(x) - \ck{x) - 



Ck{x) - [cj{x) -{y- x)Kj{y)] 



[xurjO-^i^, 



from the left 



( [Ool i-( [03i i (^ y — X 

Cn{x) - [cj{x) - {y- x)Kj{y)] 

c'"{x,y) 

c„(x) - [cj{x) - {y- x)Kj{y)] 

Cnixi) - {cj{xi) - {yi - xi)Kj{yi)} = d^^x^yi). 



(2.35) 



(2.36) 



From (2.34), (2.35) and (2.36) continuity of c" follows for any sequence {xn,yn) -* 
{x,y). We now extend the above argument to the situation when x = y which 
establishes left-continuity of c^ at {y,y). In this case we have x = C,k{y) = y- 
For this to hold we must have Ck{y) = Cj{y). Using boundedness of K^ ioi i < n 
shows that (2.34), (2.35) and (2.36) converge to each other. 

To relax (2.30) and (2.31) we successively write out Kk,Kj, . . . , until the 
assumption of the first case holds true and then, successively, apply the special 
case. 

It remains to prove continuity of C,n which we prove by contradiction. Assume 
there exist e > and y ^ such that for all 5 > there exists a. y' s {y,y + 6) 
such that \^n{y) — ^n{y')\ > e- By (2.16) the limit of ^„(?/') as y' I y exists at least 
along some subsequence and we denote it by ^n- By assumption ,^„ ^ iniy)- 

Consider first the case that ^n{y) < y and ^„ < y. Using continuity of c" we 



deduce ^"l^"i^P',f) 
Now, if 



y-£,n 



as y 



y- 



c'iLy) _^c^iUy),y) 



y-in 



y - in{y) 



(2.37) 



then we obtain a contradiction to the optimality of either ^n{y) or some Cniy') 
for y' close enough to y by continuity of c". If 



c'iLy) nUy).y) 



y-in 



y - Caiv) 



(2.38) 



we obtain a contradiction to Assumption LMii). 

We now consider the case that either ^„(y) = y or $,n = y. The case ^n(l/) < y 
and ^„ = ?/ is ruled out by condition (2.18 ) from Assumption Ani): Indeed, for the 

sequence (;.„(,) ^ ^^k^?) to be beaded we .„st Kave c"(U,),,) ^ 0. 

Recalling the left-continuity of c" at {y,y) implies Cn{y) = Ci^{y;y){y)- 
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The case ^n{y) = y and in < y follows as above by distinguishing the cases 



(2.37) and (2.38) and by recalling (2.15) and the left-continuity of c" at {y,y). □ 



Lemma 2.13 (Monotonicity of ^n)- Let n e N and let Assumption\A\hold. Th 
the mapping 



en 



in-. 



y ^ in{y) 



(2.39) 



is non- decreasing. 



Proof. The claim for n = 1,2 follows from Brown et al., 2001a. Assume by 



induction hypothesis that we have proven monotonicity of ^i, . . . , £,n-i- 



We follow closely the arguments in Brown et al., 2001a, Lemma 3.2]. Since 
in is continuous it is enough to prove monotonicity at almost every |/ ^ 0. The 
set of i/'s which are a discontinuity of ^i is a null-set, and hence we can exclude 
all such y's. In the following we fix a |/ where ^i, . . . , ^„ are continuous. 

We will first consider the case when in{y) ^ ijiy) for all j < n. By continuity 
of in it follows that there is an e > such that 

in{y) ¥= ijiy) and £ := jniv) = Jniy) for all y e (y - e, y + e) and j < n, (2.40) 

and furthermore 

iUy), y) e iUy) - e, Uy) + e) x (y - e, y + e). (2.41) 

Let /i denote a supporting tangent to c^{-,y) at inijj) which goes through the 
X-axis at y, i.e. 

li{x) = c"(e„(l/), y) + {x Uy)){D - K,{y)), 

where D lies between the left- and right-derivatives of c„ — q at in{y)- Using 
that li{y) = we can write 



D Ke{y) 



c''{Uy),y) 
y - in{y) 



Cniinjy)) - Ci{in{y)) 

y-in{y) 



Ke{y) 



and thus by (2.19) 



We also have 



D ^ 0. 



h{y + 6) = 6{D-Ke{y)). 



(2.42) 



(2.43) 



Choose S s (0,e) sufficiently small. Our goal is to prove in{y + S) ^ in{y)- 
Recall that in{y + S) is determined from y + 6 and c"'{-,y + 6) only. Since we 



13 



know that ^n{y + S) s {Cn{y) ~ e, Cn{y) + e) := / it will turn out to be enough to 
look at c"(x, y + S) only for x e (^„.(i/) — ^, ^n{y) + e)- For such an x we have 



c"(x, y + 6)^ c^ix, y)^^{y + 6-x) K,{y + 5) - {y - x) K,{y). (2.44) 
Let I2 be the supporting tangent to c^{-,y + 5) — c"-{-,y) at C,n{y), i-e. 

hix) = c-iUy),y + 5) - c'\Uy),y) + ix- Uy)){K,{y) - My + 5))- 



Hence, 



li{y + 5) + hiy + 5) 



^ 6{D-K,{y)) 

+ c"(en(Z/),Z/ + 5)-c"(a(t/),y) 

+ {y + 6 - Uy)){My) - K,{y + 6)) 



6D^0. 



(2.45) 



Now, since c"'{-,y + 6) ^c^{-,y) is linear (and therefore convex) in the domain 
J, /i + ^2 is a supporting tangent to c^{-,y + 5) at C,n{y), i-e. 



(/i + /2)(2;) < c"(x,?/ + 5) for a; e /, 

(^i + ^2)(en(y))=c"(a(2/),2/ + '5). 



(2.46) 
(2.47) 



Recall that ,^„(y + 6) is determined as the x- value where the supporting tangent 
to c^{-,y + 5) which passes the x-axis aX y + 5 touches c^{-,y + 5). Next we 
exploit the fact that ^„(|/ + (5) e J which implies that we only need to show that 
in{y + (5) ^ i^n{y) — e,^n(z/))- This follows from (2.45) which yields that any 
supporting tangent to c"(-,?/ + 6) at some ( e {^n{y) — ^,^n{y)) must be below 



the X-axis when evaluated at y + 6. We refer to [Brown et al., 20()la Fig. 7] for 
a graphical illustration of this fact. 



Now we relax the assumption (2.40). Assume that there exists a 6 > such 



that ^niv) > ^niy + S). We derive a contradiction to the special case as follows. 
Set yo := y and yn '■= y + 5 . Recall that ^„ is continuous. Now we can choose 

VO < yi < ■ ■ ■ < Vn-l < yn such that iniyo) > ^niyi) > ■■■ > ^niVn-l) > UiVn) ■ 

Set Xi := ^n{yi), i = 0, . . . ,n. Observe that by monotonicity of ^k, k < n the 
graph of ^k intersects with at most one rectangle (xj, Xj_i) x (?/i_i, yi), i = 1, . . . ,n. 
Consequently, there must exist at least one integer j such that the rectangle R : = 
{xj, Xj_i) X {yj_i,yj) is disjoint with the graph of every ^k, k < n. By construction 
and continuity of y t-^ (,n{y) R is not disjoint with the graph of ^„. Inside this 



rectangle R the conditions of the special case (2.40) are satisfied. Recalling that 
inijjj) = Xj < Xj_i = ^niyj-i) and by continuity oi y ^^ (,n{y), we can find two 
points si < S2 such that zi = C,n{si) > C,n{s2) = Z2 and {zi,si) e R, (2:2,52) e -R- 
This is a contradiction. D 
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Lemma 2.14 (ODE for Kn)- Let n s N and let Assumption \A\ hold. Then the 
mapping 

y -- K^{y) (2.48) 

is absolutely continuous and non-increasing. 



If we assume in addition that the embedding property of Theorem \2.6\ is valid 
for the first ra — 1 marginals then for almost all y ^ we have: 
IfUiv) < y then 

where K', denotes the derivative of Kj which exists for almost all y ^ and 
j = l,...,n. 

Ifin{.y) = y then 

K^iy+) = K„^f^y)iy+). (2.50) 

Proof. The proof is reported in the Appendix |A| D 

3 Proof of Main Result 

In this Section we prove the main resuh. The key step is the identification of the 



distribution of the maximum, cf. Proposition |3.4 
Let n E N. For convenience we set 



Mo:=0, Mr.= B,^, i = l,...,n, (3.1) 



where r,- is defined in Definition 2.2 



3.1 Basic Properties of the Embedding 

Our first resuh shows that there is a "strong relation" between M„ and M„. 

Lemma 3.1 (Relations Between M„ and M„). Let n s N and let Assumption\A\ 
hold. Then the following implications hold. 

Mn ^ C,n{y) =^^ Mn ^ y if Cn is strictly increasing at y. (3.2) 

For y ^ such that jn{y) j^ we have 

M,^iy) ^ Uy) > ^My)iy) =^ M.^^Uy). (3.3) 

M,4^)<1/, Mn^y =^ Mn^Uy). (3.4) 

M,^iy)^y, M,^^y^<Uy) =^ M^<Uy)- (3.5) 

If ^n is strictly increasing at y ^ and jn{y) = then the following holds. 

Mn^Uy) ^^ Mn^y. (3.6) 
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Proof. Write ] = ]n- We have 

^3(y){y) < ^n{y) ^^t{y), i = j{y) + i, • 



,n. 



In the following we are using continuity and monotonicity of ^i, . . . ,^„, cf. 
Lemma EIS and EjH 



Case j{y) ^ 0. As for imphcation (3.2) assume that ^„ is strictly increasing 
at y, Mn ^ C,n{y) and M„ < y holds. In this case M„ cannot be at the boundary 
^n- It has to be at a boundary point C,j{y') for some j < n and some y' < y. 



However, this cannot be true because ^„(y') < ^j{y') and hence case (2.7a) of the 
definition of r„ would have been triggered. 



Implication (3.3) now follows from implication (3.2) and the fact that either 



Mn = Mjiy\ (case (2.7b)) or M moves to a point at the boundary iiiy') ^ iniy) 



for some i = j{y) + 1, . . . ,n, y' ^ y (case (2.7a)) 



Implication (3.4) holds because the only possibility for M to increase its max- 



imum at time j{y), which is < y, to some y' '^ y is to hit a boundary point 
6(l/') > ^n{y) for some i = j{y) + 1, . . . , n._ 

Implication (3.5) holds because from Mj(j^) ^ y and Mj(j^) < ^„,(y) it follows 



that Mj(y) = ^i{y') < ^„(?/) < ^j{y') for some i < j{y), y' ^ y, j > j{y). From this 



it follows that M will stay where it is until time n, cf. case (2.7b). 

Case j{y) = 0. The condition M„ ^ C,n{y) implies in a similar fashion as in 
(3.2) that Mn ^ y holds. Conversely, assume that Mn ^ y holds. In this case 
Mn must be at a boundary point C,i{y') ^ Cn{y) for some i = 1, . . . ,n, y' ^ y. D 



As a application of Lemma 3.1| we obtain the following result. 

Lemma 3.2 (Contributions to the Maximum). Let n eN and let Assumption\A\ 
hold. If C,n is strictly increasing at y ^ then 

F[Mn^y]=F [Mn > Uy)] - P [M,^{y) > Uy)] + P [Mmv) > y] ■ (3.7) 

Proof. Write J = Jn- 

Case ]{y) i^ 0. Firstly, let us compute 

P[M„^l/]-P[M„^en(y)] 
^ P [M„ ^ y] - P \Mn ^ in{y).Mn ^ 1/] = P [M„ ^ t/, M„ < ^nd/)] 

= p [M„ ^ y, Mn < Uy), Mj^y) s^ y] + P [Mn ^ y, Mn < Uy), M,^y) < y] 
g F[Mn< Uy), M,(y) ^ y, Mfy) < Uy)] ■ 

Secondly, let us compute 

^[Mjiy)^y]-F[M,^y)^Uy)] 
= p [M,(,) ^ y, M,^y) ^ Uy)] + P [Mjiy) ^ y, M,^y) < Uy)] " P [Mjiy) ^ Uy)] 

1331 



p [M,(,) ^ y, M^^y) < Uy)] - P [Mn < Uy), Mj^y) ^ y, M,^y) < Uy)] 
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Comparing these two equations yields the claim. 



Case j{y) = 0. The claim follows from (3.6) and our convention cq = 0, cf. 

D 



(2.1) in Definition 2.1 



3.2 Law of the Maximum 

Our next goal is to identify the distribution of M„. We will achieve this by 
deriving an ODE for P [M„ ^ ■] using excursion theoretical results, cf. Lemma 



3.3, and link it to the ODE satisfied by Kn, cf. Lemma 2.14 



Lemma 3.3 (ODE for the Maximum). Let n e N and let Assumption \A\ hold. 
Then the mapping 

is absolutely continuous and for almost all y ^ we have: 
Ifin{y) < y then 



5P[M„ ^y]^ P [M„ ^ y] ^ P [M,„fa) ^ y] ^ 5P [M, ^ y] 



y - ^niy) 



Sy y- in{y) 

ifUiv) = y then 

P[M„>l/]=P[M,„(,)>y] 



dy 



j=3n,{y) 



(3. 



(3.9) 



Proof. Write J = jn- We exclude all y > which are a discontinuity of ^i. This 
is clearly a null-set. 



The cases n = 1,2 are true by Brown et al., 2001a . Assume by induction 



hypothesis that we have proven the claim for i = l,...,n — 1. 

If Cniy) = y then it is clear from the definition of the embedding, cf. Definition 
2.2[ that for any i < n 



Mn>y ^^ Mi>y. 

Case ]{y) i^ 0. For (5 > we have 

P [M„ ^ y + (5, Mj(j,) < y + ^^] - P [M„ ^ y, M,(^) < y\ 
= F[Mn^y + 6, M,^y) < y] - P [M„ ^ y, M,^y) < y] 
+ P [A?„, ^y + 5,y< M^^y) <y + 5\. 

' y ' 

=0 for (5 > small enough by definition of j(j/) and continuity of ^i 

For r > we define 

lj{r) := max {^^(r) : ^k{y) = ^n{y)} , 
r(r) := min {^^(r) : ^k{y) = Cn{y)} 



(3.10) 



(3.11) 
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and note that 



^Mi'^) -^ ^niy), L^ir) -> ^n{y) asr^y 



-3{y) 



by continuity of ^j at y for i = 1, . . . , n. 



(3.12) 



Let 6 > 0. We have by excursion theoretical results, cf. e.g. Rogers, 1989 



^ P [M„ ^ y + (5, Mj^y) < y] 



< 



/ r 

P[M„^|/,M,(,)<y]expl-J 



^Mi^ 



y+s dr 



-j{y)^ ' 



Now we compute for y such that ^„(y) < y 

P [M„ ^ y + (5, M,(j,) < y + (5] - P [M„ ;^ y, M,(j,) < y] 



by jIT2l P [M„ ^ y, Mj(y) < v\ 



^^"(-r-^wj-i 



and 



as5iO y-in{y) 

P [M„ ^ y + <5, M,(j,) < y + (5] - P [M„ ^ y, M,(j,) < y\ 



TTljl ( |3l3l l 



P [M„ ^ y, M,(j,) < y] 



exp (- g 



J/+5 dr 



■<--t](y){T) 



6 



by J3A2f^ P [Mn > V, M^jy) < v] 



Hence, from (3.14) and (3.15) it follows that the right-derivative of 

y^¥[Mn^y,M,<y]\^^^^^^ 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



exists. Similar arguments for 5 < show that the left- derivative exists and is 



the same as the right-derivative. Local Lipschitz continuity of (3.16) then follows 



from (3.14) and (3.15) 



Observe the obvious equality 

P [M„ ^ y] = P [M, ^y]+¥[Mn^ y, M, < y] 



(3.17) 



Taking j = j{y) in (3.17) and fixing it, we conclude by induction hypothesis that 
y >^ F [M„ > y] is locally Lipschitz continuous and hence absolutely continuous 
and its derivative reads 

dF [Mn ^ y] dF [Mj ^ y] 



dy 



dy 



+ 



j=jn{y) 



F [M,4,) ^y]-F[Mn^ y] 

y - U{y) 



Case j{y) = 0. For 5 > we have by excursion theoretical results 
F[Mn^ y + 5,Mi < y + 5\-F[Mn^ y,Mi < y] 
= F[Mn^ y + 6,Mi < y + 6] -F[Mn^ y + 6, Mi < y] 
+F[Mn> y + 6,M^ < y] -F[Mn^ y,M, < y] 



ry+o 

F [Ml E ds] 

Jv 



+P [M„ ^y,Mi< y] 



i,{s) 



exp 



y+5 



dr 



exp 



y+5 



dr 



r^^^{r) 



r-^_^{r) 
- 1 



(3.18) 



Similarly, we have 

F[Mn>y + 5,Ivh<y + 5]-F [M„ >y,Mi< y] 



r^y+S 



> 



P [Ml e dsj ^ ^ ^ exp 

■Jv 



Us) 



y+s 



dr 



+P [M„ ^y,Mi< y] 



exp 



y+5 



dr 



r-^i(r) 



r-^i(r) 
- 1 



From (3.18) and (3.19) it follows that the right-derivative of 

P [Mn ^y,Mi< y] 



y 



(3.19) 



(3.20) 



exists. Similar arguments for 5 < show that the left- derivative exists and is 



the same as the right-derivative. Local Lipschitz continuity of (3.20) then follows 



from (3.18) and (3.19). Now we can conclude from (3.17)-(3.19) applied with 
j = 1 that y »-> P [Mn ^ y] is locally Lipschitz continuous and hence absolutely 
continuous and its derivative reads 

dF [Mn > y\ 
dy 
^^dF [Ml ^ y\ dF [Ml ^ y\ {ii{y) - e„(y))+ P [M„ ^ y] - P [Mi ^ y\ 



Sy dy y- ^n{y) 

which implies by induction hypothesis 

P [M„ ^ y] dF [Mn ^ y] 



y - in{y) 



y - ^n{y) 



+ 



dy 



This finishes the proof. 



n 
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Finally, we argue that P [M„ ^ y] = Kn{y) holds for all y ^ 0. 

Proposition 3.4 (Law of the Maximum). Letn s N and let Assumption \A\ hold. 
Assume that the embedding property of Theorem 2^ is valid for the first n — 1 
marginals. Then for all y ^ we have 



F[Mn^y] =Kn{y). 



(3.21) 



Proof. The case n = 1 holds by the Azema-Yor embedding. Assume by induction 
hypothesis that 



In Lemma 



2.14 



and 



l,...,n-l. 



3.3 



we derived an ODE for Kn and P [M„ ^ ■] , respec- 
tively, in terms of i^i, . . . , Kn-i and P [Mi ^ ■],..., P [M„-i ^ ■], respectively. 
These ODEs are valid for a.e. y ^ 0. By induction hypothesis both drivers of 
these ODEs coincide everywhere and hence the claim follows from the boundary 
conditions 



Kn{y) ^0 as y 

P [M„, ^ y] ^ as y 



00, 
00, 



Kn{y) ^1 as y ^ 0, 

P [M„ ^ y] ^ 1 as y -> 0, 



absolute continuity of Kn and P [M„ ^ •] and the fact that the ODE 

¥[Mn^y]-Kn{y) 



¥[Mn^y]-Kn{y) 



has unique solution given by 0. 



y - in{y) 



P [M„ ^ 0] - ir„(0) = 0, 

n 



3.3 Embedding Property 



, r„ from Definition 2.2 



In this subsection we prove that the stopping times ri, 

embed the laws /^i, . . . ,/i„ if Assumption [A] is in place. More precisely, given 



Proposition |3.4| above and by inductive reasoning, to complete the proof of The- 
orem 



2.6 we only need to show the following: 



Proposition 3.5 (Embedding). In the setup of Theorem 2.6 we have 
and {Br„/,t)t>o ^■^ '^ uniformly integrable martingale. 



(3.22) 



Proof. The case n = 1 is just the Azema-Yor embedding. By induction hypoth- 
esis, assume that the claim holds for all i ^ n — 1. 
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We claim that ^„ ranges continuously over the full support of /i„. This is be- 
cause, firstly, we know from Lemma 2.13| that C,2, ■ ■ ■ ,^n are continuous. Secondly, 
we have by using l^^ < l^- that 

■ , c"(C,o) ^. , . ) c„(C)-q(C) , ,^,,, ( c„(/^„ 
ml — — ^ mt mm < : , h Kj(0) > a 



C^O — C C=£0 1s;i<ri 



0-C 



O-L 



^0 



which shows that ^n(O) = /^„. Furthermore, by using r^^ ^ r^^ we have from 
^J^ and (|2l6l) that 



^.nVflnj 



fin- 



Let y > be such that ^„ is differentiable and strictly increasing at y, ^n{y) is 
not an atom of neither /i„ nor Hj„(y) and y is not a discontinuity of ^i. Note that 
for such a y equation (A. 6) holds because of (A. 7). Applying previous results we 
obtain 



p [Mn ^ Uy)] - P [M^y) ^ Uy)] + P [M^y) ^ y] 



Lemina 13.21 
Prop. [3]4] 



P [M„ ^ y] 

Kn{y) 

- C'n{in{y)) + Cj„(2;)(^n(2/)) + ^j„(2/)(2/)' 

which implies by induction hypothesis that 

P[M„ ^ en(l/)] = ~(^rXUy)) = /in(Kn(y),<X.)). 

We have matched the distribution of M„ to /x„ at almost all points inside the 
support. The embedding property follows. 



Now we prove uniform integrability by applying a result from Azema et al.. 



1980 which states that if 



lim xP[|5|^^ ^x] =0 



(3.23) 



then (-Br„At)i^o ^^ uniformly integrable. 

Let us verify (3.23). Set H^ = inf {t > : i?t = x]. We have (here ^^^ denotes 
the left-continuous inverse of ^j) 



f\\B\ ^xl ^P 



H^,. < H 



maxis;„5j ^{-^) 



+ P [S,„ ^ x] 



maxi^„ ^, 



-1/ 



X + maxi^n^i (-x) 



+ KJx). 
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From the definition of ^„, cf. (2.4), and tlie properties of bi, cf. (2.11 ) we have 



niax^^ '^{—x) ^ niax6j(— x) — > 



i^n 



i^n 



and hence, recalhng the definition of /x^^ in (2.23) 



c (b~^(x)) 
lim xF [ISL ^ x] < lim xKJx) < hm x " , )/[ = lim a;/x5?^([a;, oo)) = 0. 



»oo X — b~^{x) X' 



This finishes the proof. 



n 



4 Discussion of Assumption |A| and Extensions 

In this section we present a simple example of probability measures /ii,/i2,/^3 
which violate Assumption lAFii) and where the stopping boundaries ^1,^2,^3, ob- 



tained via (2.4), fail to embed (/Xi,/i2,/i3). 

Motivated by this counterexample we sketch how an iterated Azema-Yor type 
embedding in the three marginal case n = 3 looks like, i.e. in the absence of As- 
sumption |A](ii). This will lead us to redefine ^3 in certain regions. More precisely, 
we show how to modify the optimization problem from which ^3 is determined 
in order to obtain the embedding property. We will also point out some ma- 
jor differences of this (general) embedding to the embedding in the presence of 
Assumption lAFii) . 



4.1 Counterexample for Assumption A(ii) 

In Figure |4.1| we define measures via their potentials 



Ufi: 



Ufi{x) :-- 



\u — x\ fi{du). 



We refer to Obloj, 2004, Proposition 2.3] for some properties of Ufi. 
The measures are given as 



/^l(-l) = 


2 
^3' 


/^i(2) = I, 






/^2(-3) = 


2 

~~7' 


fl\ 18 

/^2 7^ = ^' 

\2J 35 


/i2(3) = 


1 

" 5 


/"3(-3) = 


2 

~~7' 


/^3(-2) ^ -, 


/i3(3) = 


16 
" 35 



(4.1^ 



(4.2) 
(4.3) 
(4.4) 



Observe that the embedding for (/ii,/i2,/i3) is unique: We write Ha^ for the 
exit time of [a,b] and denote Ha^t ° (^t '■= inf {t > t : Bt ^ {a,b)}. Then the 
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Figure 4.1: Potentials of /Ui, /i2, /^3, t^2 and z/3. 



embedding (ri, r2, Ts) can be written as 



Ti = i^-1,2, T-2 = i:f_3 1 o 6*^^ 



+ H1.0 9^^ 



, Tg = i:f_2,3 O 61^ 



(4.5) 



As mentioned earlier, our construction yields the same first two stopping 



boundaries as the method of Brown et al., 2001a . In this case, cf. Figure 



4.2 



-1 if ye [0,2), 
y else. 



'-3 a ye [0,1 



Uy) ■■- 



ifl/6[i,3) 



2^' 

2)3), 



y else. 



This already shows that our embedding fails to embed fi2'- while the embedding 
for fii is valid by the Azema-Yor type embedding, the embedding for /i2 is invalid. 
To see this one just has to compare the stopping boundary ^2 in the Definition 



of T2 with (4.5). In Section (4.2) we will recall from Brown et al., 2001a how 



the stopping time T2 has to be modified. 

However, the point we want to make, is that the embedding for fi^ fails 



and cannot be achieved by using ^3: the optimization problem (2.4) does not 



23 




Figure 4.2: We illustrate the (unique) boundaries r]i,r]2,i]3 required for the em- 



bedding of (/ii,/i2,/i3) from (4.2)-(4.4) and the stopping boundary ^3 obtained 
from ( 2.4[ ). We already know that 771 = ^1 and t]2 = ^2- In order to ensure the em- 
bedding for fi2, the mass in the region (—1) x (^, 2) is redistributed to (—3) x (^, 2) 
and (0.5) x (^,2). Note that the case ^2(2/) = y = ^2(1/) is possible and required 
to define the embedding. After the embedding for (/ii,/i2), the Brown et al., 



2001a embedding, our goal is to define the embedding for fi^ on top of that "in 



the same fashion". In this example 773 is implied directly by (4.5). In Section 4.2 
we show how to obtain it in general. 



return the third (unique) stopping boundary which is required for the embed- 
ding of (/ii,/i2,/i3)- Indeed, for sufficiently small y > ^, in the region ( < 
min(^i(|/),^2(?/)) = "1 "we are looking at the minimization of (^ >— > ^^^^ which is 
attained by ^3(1/) = — 3 < — 2 since fi^ has an atom at —3. Consequently, there 
will be a positive probability to hit —3 after T2. This contradicts (4.4). 



This example does not contradict our main result because AssumptionL4[ii)(a 



is not satisfied for i = 2 and y 
but C2 (^) = ci (^) holds. 



2' 



where C = —3 minimizes the objective function 



4.2 Sketch for General Embedding in the Case n = 3 



In the example of the measures (/xi, /i2, /is) from (4.2)-(4.4), we note however that 
the (unique) embedding resembles something which might be termed "iterated 
Azema-Yor type embedding" although this embedding does not satisfy the rela- 
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tions from Lemma 3.1 Consequently, one might conjecture that a modification 



of the optimization problem (2.4) and a relaxation of Lemma 3.1 might fix this 



problem. We now explain in which sense this is true. The purpose of this section 
is to outline the main arguments and hence we sacrifice occasionally on rigor. 
In order to understand the problem in more detail, we need to recall from 



Brown et al., 2001a how the embedding for fi2 looks like in general. It reads 



BHR 

^2 



in2'^UBr 



and6(5.J<6(^. 



T2 else. 



(4.6) 



where T2 is some stopping time embedding whose existence is established in 
Brown et al., 2001a and corresponds to the case where /ii and fi2 have in- 
tervals which have same mass, mean and are in convex order. One important 
property of T2 is that it does not change the current maximum, i.e. -B,-^ = -B,-' . 
In general there will be infinitely many such stopping times Xg. Although this 
is not true for (/ii,/i2,yU3) because their embedding is unique, we illustrate this 
claim by the measures (yUi, z/2, z/3) which are defined via their potentials in Figure 



4.1 This example also suggests that the atomic nature of the measures (which 
results in discontinuous stopping boundaries) is not the main problem but rather 
the intervals with same mass, mean and convex order. 



Let ^1 and ^2 be defined as in (2.4) but now we let T2 : 



et al., 2001a embedding, cf. (4.6). Also set M2 = B. 



^BHR ]^g |-]^g Brown 



■'T2- 



Now our goal is to define an embedding t^ for the third marginal on top of 

We still want 



the [Brown et al., 2001a embedding in a situation as in (/ii, 1/2, z/3 

to define our iterated Azema-Yor type embedding through a stopping rule based 

on some stopping boundary ^3 as a first exit time. 



inf 



^3 : = 



t ^ r^^HR 



T, 



BHR 



■.Bt^UBt)} ifi?.BHR>e3(A 

else, 



-BHR 

2 ' 



(4.7) 



and prove that this is a valid embedding of ^^. We observe that now the choice of 
Tg in the definition of T2 = r^^^ matters both for the embedding and the distri- 



bution of the maximum B^^. Similarly as in Brown et al., 2001a we expect that 



this will be only possible if the procedure which produces ^3 yields a continuous 
^3. We restrict to this case because otherwise an additional step is likely to be 
required. 



With this, a more canonical approach in the context of Lemma 3.1 is to write 



P[M3^1/] =F M^^Uy) 



+ "error-term". 



which we formalize in (4.19). As it will turn out, this "error-term" provides a 
suitable "book-keeping procedure" to keep track of the masses in the embedding. 
Then we proceed along the lines of the proof of our main result. 
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As mentioned above, we do not present a rigorous proof for the embedding 
property but rather outhne the main ideas. For simphcity, we also assume that 
^2 has only one discontinuity, i.e. S_ := ^2(1/—) < ^2(1/+) '■= 2+ for some y ^ 0. 
Denote y := ^r^(H+). 



4.2.1 Redefining ^3 and K^ 

Define the following auxiliary terms, 

f^''(C,y;rD:=E[F(C,l/;r^)]. 
Note that for y ^ y, 



dC 



-(C,i/;r^)=P[Mie[2/,y],M2<C] 



and 



df 



iAY 



dy 



-iC,y;^2] 



-E 



^{Miedy,M2<C} 



dy 



C + E 



-{A/iedy,A^2<C} 
d^ 



Mo 



C-a(C,i/;^2)) 



X P [Ml E dy, M2 < C] 



dy 



where 



a 



KC,l/;r^):=E[M2|Mi=l/,M2<C], 
/3(C,l/;r^):=E[M2|Mi=y,M2^C]- 

With these definitions we have by the properties of T2, 



<iiy) 



P [Ml e di/] 
dy 



dy 



dy 



(4.8) 
(4.9) 

(4.10) 



(4.11) 



(4.12) 
(4.13) 



(4.14) 



We now redefine ^3 and K3 from (2.4) and (2.6), respectively, and denote the 



new definition by ,^3 and K3. To this end, introduce the function 

C3(C) - f'-^iC, y, r',) if f ^Y(C, y; r',) > and C < S+, (4.15a) 



^ ^ ' c^(C,l/) else. (4.15b) 



Now set 



6(2/) := argmm — 



(4.16) 
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and 



KM :-- 



c^{^3{y),y) 



y - 6(2/) 

The first order condition for optiniality of ^3(1/) reads 



(4.17) 



iAY 



Ks{y) + c'^iisiy)) ^(^(z/),!/; r^) = 0. 
We stress that now the quantities ^3 and K^ depend on Xg. 



(4.18) 



4.2.2 Embedding Property 

In the following we assume that ( and y are such that /'^"^(C?!/; ^2) > and 
C < S_|_. If this is not true the arguments from Sections [2] and [3] apply. 



By definition of the embedding in (4.7) we have, assuming ^3 is non-decreasing, 

P [M3 ^ y] = P 



P 
P 



M3 ^ Uy) 
Ms > Uv) 

Ms ^ ls{y) 
Ms ^ ls{y) 
We observe the obvious equation 



+ P 
+ P 
+ P 
5/ 



+ 



5C 



Ms^y,Ms<is{y) 
Mi^y,M2<Uy) 
M^e[y,y],M2<Uy) 

iAY 

fe(y),i/;7-2). 



(4.19) 



P [M3 ^ 2/] = F [Ml < y, M3 ^ y] + P [Ml 6 [y, y), M3 ^ y] + P [Ml ^ y] 

(4.20) 



and note that 

aP[Mie[y,y),M3^m] 



P 



dy 
M2>Uy),M,Edy 



p{i3{y),y;^2) ■-- 



m=y 



dy 



P 



Hy o e,, < if^-3(^) o ^,, |M2 > fsd/), Ml = y 



(4.21) 



where as above H^ := inf [t 1^ : Bt = x) . 
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Hence, by (|420 



, (4.21) and similar excursion theoretical results as in the 



proof of Lemma 3.3 

d r - T FlMs^y] ~ ^ ^ ,, 5P [Ml ^ yl P [Ml ^ yl 

-p M3 ^y\ = — ^ I /^ + p{Uy). y. r + ^ , ^^ + ^ - , ^^ 
^2/ y-^sd/) ^'^ y-Uy) 

F[M,^y\ ~ ^,{y) - ^,{y) dF [M, ^ y] 

+ p{^3{y),y;r^) 



y - 6(y) 



y - 6(y) 



dy 



(4.22) 



KiyJ 



Hence, by similar calculations as in (A. 17) 

Ksiy) '-^iUy),y;r^) 



y - Csiy) 



y - ^siv) 



jint Ksjy) ^ 6(y) -a{^3{y),y) 



p 



Misdy,M2<Uy) 



y - 6(y) 



y - i?>{y) 



dy 



Hi ^3(1/) , Uv)-iM^{^^^M I /3(i/)-e3(y) 



p 



Miedi/,M2^e3(l/) 



y-iz{y) y-iM ^y y^Uy) 

jm KM Uy) Uy) ^^ [^i > y\ 



y - i-Ay) y - iz{y) 



dy 



+ p(6(2/),2/;7-2) 



Consequently, as in and together with Proposition |3.5[ for all y, 

kM=^{M^^y\ 



and 



-4(6(2/)) = P M3 ^ Uy) 



dy 

(4.23) 

(4.24) 
(4.25) 



If ^3 happens to be discontinuous, then we cannot guarantee the embedding 



property for jj^. However, we expect that similar arguments as those in Brown 



et al., 2001a, Section 3.5] will yield a stopping time Xg which does the job for the 
jump-intervals of ^3. We do not pursue these matters here. 



A Appendix: Proof of Lemma 2.14 



In order to prove Lemma 2.14 we require to prove, inductively, several auxiliary 
results along the way. We now state and prove a Lemma which contains the 
statement of Lemma 12.141 
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Lemma A.l. Let n e N and let Assumptions^ hold. Then the mapping 

y - KJy) (A.l) 

is absolutely continuous and non-increasing. 



If we assume in addition that the embedding property of Theorem \2.6\ is valid 
for the first n — 1 marginals then for almost all y ^ we have: 
Ifin{.y) < y then 



Kiy)^,^^^-K,y,iy)^''-'^'^'^ 



(A.2) 



y-Uy) '"'''"' y-Uy) 

where K', denotes the derivative of Kj which exists for almost all y ^ and 
j = l,...,n. 

Ifiniy) = y then 

K^{y+) = Kj^^y){y+). (A.3) 

For a; > the mapping 

c":(x,oo)^M, y^c'^^x.y) (A.4) 

is locally Lipschitz continuous, non- decreasing and for almost all y > 

do'' 



dy 



{x,y) 



^■=6i(s/) 



Kjr.{y){y) + {y- ^n{y))K'^„(y){y)- 



(A.5) 



// the mapping c'^{-,y) is differentiable at S,n{y) o-nd ^^{y) > then for almost 
ally^O 



Kn{y) + c'^iUy)) " cAUy)) - K,{y) = o 



(A.6) 



for j = jn{y) and j such that n > j > jn{y) and C,n{y) = C,j{y)- In the case of 
non-smoothness ofc^{-,y) at C,n{y) we have 



CXy) = 



(A.7) 



fory such that the slope of the supporting tangent to c'^{-,y) at^n{y) which crosses 
the X-axis at y does not equal the right-derivative of c'^{-,y) at S,n{y)- 

Proof. We prove the claim by induction over n. The induction basis n = 1 holds 
by definition and Brown et al., 2001a, Lemma 2.6]. 

Now assume that the claim holds for alH = 1, . . . , n — 1. 

Induction step for c". We have 



c"(x,|/ + 6) - c"{x,y) = - [c^„(x;y+s){x) - {y + S - x)K,^(^x.y+s){y + S)] 

+ [ct„{x;y){x) -{y- x)K,„^x.y){y)] . 

29 



(A.8) 



Firstly, consider the case when there exists a S' > such that for all |(5| < (5' 



we have tn{x; y) = tn{x', y + S). Equation (A. 8) simplifies and we have 

c"(a;, y + S)- c"(x, y) = {y + 6 - x)K^^(^^.y) {y + 6) - {y - x)K^^^^.y) (y) 

= {y + S- ^tn{x;y){y))K^„^x;y){y + S) - {y - ^t„(x;y){y))K,^ix;y){y) 
+ {X - ^^n{x■,y){y)) [K^„{x■,y){y) " K,^{x;y){y + S)] 

^ c-(^^^)fe„(.;,)(i/), y + 6)^ c^"^''^''\^^Ax■,y)iy), y) 

+ {y- itr^{x;y){y)) [K,^(x;y){y) " K,^{x;y){y + S)] 

< max \d{Uy),y + s) - c\Uy),y) + iy- Uy)) [K^{y) - K,{y + s)] } 

i<n I J 

< const(|/) • \6\ (A.9) 

by induction hypothesis and where const (?/) denotes a constant depending on y. 
A similar computation shows for \6\ small enough 

d'{x,y)-c''{x,y + 6) 

+ {x - ^t„(x;y){y + S)) [K^^(^^,y){y + 6)- K^„(^^,y){y)] 

const(|/) • \6\ if (5 < 0, 
^0 if (5 ^ 0. 



(A.IO) 



Monotonicity of c"(x, ■) follows. Equation (A.IO) together with (A.9) imply the 



local Lipschitz continuity. Plugging x = ^niy) into (A. 8), direct computation 



shows that (A. 5) holds 



Secondly, consider the case when z„(x; ■) jumps at y. Note that this is only 
possible when x satisfies 



X = ii^{x;y-5){y) for (5 > small enough. 



(A.ll) 



i.e. when x = ^k{y) for the index k = tn{x; y — 6) > jn{y)- By (2.9) there exists a 
6' > such that inix; y + 6) = in{x] y) for all < (5 < 5' . Hence, for 5 > small 
enough, \c^[x,y + 5) — c"(x,|/)| has the same upper bound as in the first case. 
Monotonicity of c"'(x, ■) follows. 



Furthermore, for 5 > we have for the x in ( A.ll ) that «„(x; y — 5) > tnix; y) 



holds. For notational simplicity we only consider the case 

tt„{x-y-5){x;y - 5) = in{x;y). 



The general case follows by the same arguments. We deduce from (A. 8) the 
following two equations, 

d'{x,y- 6) - d'ix^yf^iy -6- x)K,^^^^,^^^^^(^^,y^s){y - 6) - {y - x)K,^i^^.y){y) 

= {y^S- x)K,^^^.y){y -5)~{y~ x)K,^^^.y){y) (A. 12) 
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and 



c"(a;,y)-c"(x,y-(5) 



(|AlT|l 



-{y-S- x)K,„^^.y_s){y -S) + {y- x)K,^^^.y_s){y)- 



Now the local Lipschitz continuity of c"(x, •) follow from the above two equations 
and the first case. 



We prove (A. 5) by computing the required right- and left-derivative of c"(x, 
at X = $,n{y)- The right-derivative is simply, using (2.9[) and (A. 8), 



and the left-derivative is, writing k = in{in{y)\y~) > in{^n{y);y) = Jniy), 
1 



lim ^ ( - Cki^niy)) + {y + S- in{y))Kk{y + 5) 

(5tO \ 

+ Cj4y)(^n(y)) - (y - in{.y))Kj^(^y){y) 

«"(^)=«'=(^) lim \{{y + 5- Uy))My + s)-{y- Uy))KM 

(5tO V 



Ku{y) + {y~ Uy))K'M ^ K.^ivM + {y~ UyM^ivM (A-14) 



by induction hypothesis. So the two coincide for almost all y > 0. 

Induction step for Kn- A straightforward computation shows that the map- 
ping y 



^^'"^^ is non-increasing and hence for 5 > 



Kn[y + 0) = mf — < mf — < mf — = Kn{y) 

c.^y+5 y + b-C C^y y + d-C C^yy-C 

proving that Kn is non-increasing. 

Using that c"(x, ■) is non-decreasing and that ^„ is continuous, local Lipschitz 
continuity of Kn now follows from 



KM ^ £!M^±^ ^ c"(Uy + S),y^S) ^ ^^^^ ^ _, 



y - in{y + 5) 



y - iniy + S) 



6 



y - in{y + 5) 



if in{y) < y and if ^n{y) = y, recalling (2.15), we have 

'c„(C) -Q„(C;y+)(C) 



Kn{y+] 



inf 



Lemma [2713 



inf 

y^«y+) 



{y+) C 

C«(C) - C^^{C■,y+){0 

{y+) C 



+ K^AC;y+)iy) 

+ ^*n(C;2;+)(l/+) [ = K^„{y■,y)iy+) 



Kin{y)iy+), 
and local Lipschitz continuity of Kn follows by induction hypothesis. Equation 



(A. 3) is also proven. 
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Local Lipschitz continuity of c"( 
(2.8), the fact that the functions Cj 



y) follows from the properties of «„, cf. 
i = l,...,n, are locally Lipschitz and a 
similar expansion of terms in the case when ^n{y) = ^iiv) for some i < n Let 
k = '>'ni^n{y) + ', y)- This also shows that for any y ^ kinks of c"(-, y) come from 
kinks in C4,i = 1, ... ,n. 



In order to prove (A. 6) we first exclude all y ^ such that ^niv) is an 
atom of Cn and Ck and Cniv) > 0- Amongst all y s {^^ > 0} this is a null- 
set. By assumption c'^{-,y) is different iable at ^n{y)- Recalling the equations 
(2.34)-(2.36), a direct computation proves (A. 6) for k = in{^n{y) + ',y)- Now we 



want to apply the induction hypothesis to c . By choice of y we have that Ck is 
differentiable at ^n{y) = (,k{y), i-e. fik does not have an atom at ^n{y)- Hence, 
by the assumption that the embedding for the first n — 1 marginals is valid we 
cannot have Ckiy) = (except on a null-set because otherwise the embedding 
would fail). By (A. 7), c^{-,y) therefore has to be differentiable at C,n{y) = ^k{y)- 
This shows that we can indeed inductively apply (A. 6) to deduce for j = jn{y) 
and j such that n > j > jn{y) and C,n{y) 



Q{y) the following equation, 
= K^iy) + c'^iUy)) - c'.iUy)) + Kkiy) 

= K^{y) + d^{Uy)) - c'kiUy)) + c'kiUy)) - c^iUy)) + K,{y) 
= K^y) + c'^iUy)) - c'.iUy)) + KM- 



Equation (A.6) is proven. 



For later use we note the equation 

Cniiniy)) - Cj„iy)i^niy)) 



y - ^n{y) 



+ Kj^(y){y) = Kn{y) 



Cniinjy)) -CkjCnjy)) 

y - ^n{y) 



+ Kk{y) 
(A.15) 



for k such that n > k > jn{y) and ^„(y) = ^k{y)- 

Finally, we prove the claimed ODE for Kn in the case C,n{y) < y- For almost 
all 2/ ^ we have 



Kiy) = lim-T 

5^0 

= lim - 

5^0 6 



c^iUy + s),y + s) c^{Uy),y) 



y + S- in{y + 5) 
1 



y - in{y) 
1 



y + ^- ^n{y + S) y- in{y) 

c"(en(2/ + 5),l/ + <5)~c"(a(l/),2/) 



c"(a(i/ + (5),y + (5) 



i'n{y) 1 
y - in{y) 



Kn{y) + 



y - in{y) 
1 



y-^n{y) v^o 



^■^ c"(a(y + ^),y + ^)-c"(en(l/),y) 



The main technical difficulty comes from the possibility that ^n{y) = ^k{y) for 
some k < n. We present the arguments for this case and leave the other (much 
easier) case, to the reader. 



32 



By assumption the last limit exists and hence we can compute it using some 
"convenient" sequence Sm i where 6m is such that jn{y + Sm) = I for all m e N. 
Recall ( 2.9[ ). For 5m small enough such that ini^niv)', y + 5m) = Jn{y) "we obtain 



C"(^n(y + 5m), y + 5m) " C^i^niv), V + 5m) 
-- Cn{in{y + 5m)) ^ Cl{^n{y + 5m)) + {y + 5m ^ Ca{y + 5m))Ki{y + 5„ 
-Cni^niv)) + Cj„{y){^n{y)) - {V + 5m - inijj)) K j^^y){y + 5m) 



l |Xl5l l 



Cn{in{y + 5m)) ^ Cl{^n{y + 5m)) + (^ + <^m - ^n(Z/ + 5m))Ki{y + 5„ 



K 



3n{y) 



(y)) 



-CniUy)) + ciiUv)) -{y~ Uy)){Ki{y) 

-{y + 5m - in{y))Kj^^y){y + 5m)- 

From this we obtain for almost all y ^ by using the induction hypothesis 

jj^ C^'i^niy + 5m),y + 5m) - C"(^n(y), y + 5m) 
m—i-co 5m 

= ay+)[cniUy)+) c[{Uy)+) Uy) 

+Ki{y) + (y - Uy))K'iiy) - K,^(y)iy) - (y - Uy))Kiy)(y) 
g ~ay+)Kn{y). (A.16) 

1A.21 



Together with (A. 5) this yields in the case when c"'{-,y) is differentiable at 

Uy) 



Kiy) 



ay) 



y - in{y) 
Kn{y) 



1 / dd^ 

Kn{y) + y_^^^y) \^-K^{y)en{y) + j^iUy),y) 



+ 



y - ^y) y - ^y) 



K 



3n{y) 



{y) + {y-Uy))K'Uy))- (A.17) 



In order to finish the proof we just have to establish that (A. 17) also holds in the 
case when c"'(-,|/) is not differentiable at C,n{y)- 



To this end, we first argue that (A.16), and hence (A. 17), remains true in 



the case when c"(-,y) is not differentiable at iniy), but when the slope of the 
supporting tangent to c"{-,y) at C,n{y) which passes the x-axis at y equals the 
right-derivative of c"'{-,y) at C,n{y)- In that case, denoting k = in{^n{y) + ] y), 

c'n{Uy)+) c',{Uy)+) - Kiy) = ~~Kn{y)- (A.is) 

Recall the sequence {5m)- We do not necessarily have k = jn{y + 5m) = l- 
Nevertheless, we argue that 

Cn{Uy)+) cAUy)+) Uy) = ^km (a.i9) 

holds. We can safely assume that ^^(i/+) > (in the other case the conclusion of 



(A. 17) remains true). Also it is enough to consider the case when k > jn{y + 5) for 
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all (5 > sufficiently small (otherwise we may consider an alternative sequence 
{6m) where k = jn{y + <^m) for all m and (A. 19) would follow from (A.18)). 
Consequently, Ckiy~^) ^ C(l/+) > 0- Then, since by induction hypothesis (A. 2) 
holds true for k, we must have, for almost all y that 

CkiUy)+) criUy)+) - KAv) = -Kk{y) (A.20) 

holds for some j ^ /. We can iterate this argument until j = I. Then, combining 



these results we conclude by (A.18) and (A.20) that indeed (A. 19) holds 



Now we consider the case of non-smoothness of c"'{-,y) at ^niy) and where y 
is such that the slope of the supporting tangent to c^{-,y) at ^n{y) which crosses 
the X-axis at y does not equal the right-derivative of c^{-,y) at C,n{y)- We show 
that in this case we have for sufficiently small 6 > 0, 

Uy) = Uy + S) and hence C(y+) = 0, (A.21) 



which implies that (A. 17) holds as well. 



To achieve this we place a suitable tangent to c"(-,?/) at iniy)- Since, by 
assumption, c^{-,y) has a kink at ^n{y) we have some flexibility to do that. 



Recalling the tangent interpretation of (2.14) we know by choice of S,n{.y) that 
we can place a supporting tangent to c'^{-,y) at ^„(y) which passes through the 
X-axis at y. Alternatively, by choice of y, we can place a tangent to c"'{-,y) at 
C,n{y) which crosses the x-axis at some y + 6 > y. This implies that 



argmm — '-"" = ^„(i/). 



(A.22) 



Assume flrst that C,n{y) < y- Denote k = in{(,n{y) + ] y)- For simphcity of the 
argument let us also assume that ^n{y) = C,k{y) ^ Ci{y) ^^i all i i^ k,n. Also 
denote j = ]n{y) = in i^njy) ; y)- 

Now we will use ( |A.22 ) to deduce (A.21). By continuity and monotonicity 
of ^n we have for 5 > small enough that ^nijj) ^ inijj + 6) < ^n{y) + ^ < y 
for some e = e{6) > 0. By taking 6 small enough we can also assume that 
k = max^^^^(j,)+e ^n(C; y + ^)- Then we have 

■ , c^iC,y + S) ^ ., c^iCy) , ., c^iC,y + 6)-c^{C,y) 
ml — ^ ml + ml 

C^y+S y + d - C ^n{yHC<n{y)+e y + - ( ?n(y)=£C=£Cn(»/)+e y + - ( 



(A.23) 



As for the flrst inflmum in (A.23) we know from (A.22) that it is attained at 
C = ^n{y)- Now we will show that the second inflmum in (A.23) is also attained 
a-t C = ^n{y)- To this end consider the following estimate for ( > C,n{y), 

= -c*„(C;j/+5)(C) + c^„(C;^/)(C) - (z/ - OK^„iC■y){y) + {y + S - ()K,^{c,y+6){y + 5) 

^-{y - C)^v(c;«)(C;s/)(l/) + {y + S~ C)^«„(C;?/+<5)(l/ + ^) 

^-{y - C)K,{y) + {y + 6- C)K,{y + 5) =: /(C, y; 6). (A.24) 
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Since l{-,y; S) is non-decreasing and non-negative, we deduce that 

argmm — — = ^„(y). 



Finally, because at C = iniy) there is equality in (A. 24) we can conclude 



argmm — — = e.n[y) 

c^y+6 y + d-C 

as required. 

In the case when ^n{y) = y we obtain by (2.9) for 5 > sufficiently small that 
Zn{y; y + 5) = in{y] y) and hence 



arg mm — — = arg mm 

C<y+5 y + - Q y^(<y+5 



( \ 

——1. 7 + ^*n(C;y+5)(Z/ + '^) 

Z/ + ^~s ^ V ' 

The proof is complete. D 
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